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Abstract 



We introduce a new model of a stochastic sandpile on a graph G containing a sink. When 
unstable, a site sends one grain to each of its neighbours independently with probability p € 
(0,1]. For p = 1, this coincides with the standard Abelian sandpile model. In general, for 
p E (0, 1), the set of recurrent configurations of this sandpile model is different from that of the 
Abelian sandpile model. We give a characterisation of this set in terms of orientations of the 
graph G. We also define the lacking polynomial Lq as the generating function counting this set 
according to the number of grains, and show that this polynomial satisfies a recurrence which 
resembles that of the Tutte polynomial. 

1 Introduction 

In this paper, we analyse a stochastic generalisation of the Abelian sandpile model (ASM). Infor- 
mally (we provide a formal definition later), the ASM operates on a graph where each vertex has a 
number of 'grains of sand' on it. At every unit of time, another grain of sand is added at a random 
vertex v. If this causes the number of grains at v to exceed its degree, v topples, sending one 
grain to each of its neighbours. This may cause other vertices to topple, and we continue until the 
configuration is stable, i.e. no vertex can topple anymore. A special vertex, the sink, can absorb 
any number of grains and never topples. It is possible to show that eventually this model will be 
trapped in a set of configurations, called recurrent configurations. 

This model arose from work by Bak, Tang and Wiesenfeld [21 [3] and was named and formalised 
by Dhar [H] . It displays a phenomenon known as self-organised criticality |16[ [S] , in which char- 
acteristic length or time scales break down in the 'critical' steady state. When this happens, the 
correlation between the number of grains at two vertices obeys a power-law decay, as opposed to 
an exponential decay often found in models away from criticality. Likewise, the average number of 
topplings that result from a single grain addition also obeys a power-law distribution. In this sense, 
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the model is 'non-local', as grains added at a vertex may have an effect on vertices that are far 
away. Physically, this model (and self-organised criticality in general) has been used in applications 
as wide as forest fires [E], earthquakes PQ and sediment deposits [H]. 

Mathematically, the ASM has been heavily studied, and we shall not list out all the references 
here. We refer interested readers to Dhar's papers [61 [7] and to the excellent review on the subject 
by Redig [12]. Some relevant results to our current work discussed in [SIE] include: 

• The number of recurrent configurations is equal to the number of spanning trees of the graph. 

• There is an algorithm (called the burning algorithm) which determines if a given configuration 
is recurrent or not. This finds, or establishes the non-existence of, a subgraph not including 
the sink on which the configuration is stable. This algorithm constructively establishes a 
bijection between recurrent configurations and spanning trees. 

• In the steady state of the model, each recurrent configuration is equally likely. 

In the ASM, the only randomness occurs in the vertices that we add grains to. We introduce a 
variation on this model, where the topplings themselves are also random. More precisely, we fix a 
probability p G (0, 1], and when a site is unstable, each neighbour independently has a probability 
p of receiving a grain from the unstable site. In this way, an unstable site may remain unstable 
after toppling but, as in the original model, the process continues until the configuration is stable. 
If p = 1, this is identical to the ASM. 

Although this new model appears similar to the ASM, a closer inspection reveals some quali- 
tative differences. In particular, the model will again become trapped in a set of recurrent config- 
urations, but this set is not equal to the set of recurrent configurations in the ASM. Furthermore, 
each recurrent configuration is not equally likely, and the steady state measure now depends on p. 
The aim of this paper is to study the behaviour of this new model, particularly in these respects. 
We prove a characterisation of the recurrent configurations in terms of orientations of the graph 
edges, and provide a "Tutte polynomial-like" formula which counts these configurations in terms 
of their numbers of grains of sand. 

The stochastic sandpile model (SSM) we propose is an appropriate generalisation of the ASM, 
due to the aforementioned relation between the Tutte polynomial and a counting of the recurrent 
configurations. A similar result has been proved by Lopez [TU] for the ASM; in this model, the 
lacking polynomial we define later has an expression using the "standard" Tutte polynomiaQ See 
also Cori and Le Borgne [5] and Bernardi [5] for combinatorial explanations of this fact. Although 
several instances of random sandpile models have already been introduced in the literature (see 
for example Dhar Section 6], Kloster et al. [SJ, Manna as far as we are aware none of 

1 Indeed, one can pass from their "level polynomial" to the lacking polynomial by a change of variables and the 
multiplication by a monomial factor, since the level of a configuration is r\ v — \E\ if the number of lacking particles 
is 2\E\ - l-J2Vv- 
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these models have any known link with Tutte-like polynomials. Variations on the basic sandpile 
model also include an asymmetric version [15] , which devised an analogous algorithm to the burning 
algorithm called the script algorithm. 

In Section [2j we formally define the ASM and introduce our stochastic generalisation of it (the 
SSM). In Section [3j we present our results on the SSM, giving detailed proofs in Section |4j Finally, 
we offer a brief conclusion in Section [5j 

2 The model 

In this section, we formally define our model and its associated notation. We start with a brief 
review of the established Abelian sandpile model. 

2.1 The ASM 

We first define the class of graphs G = (V U {s},E) that underly the model. G must be finite, 
unoriented, connected and loop-free. It may have multiple edges, and it contains a distinguished 
vertex s that we call the sink. The set of these graphs is denoted Q. We use the notation u ~ v to 
denote that u and v are adjacent in G, i.e. {u, v} G E. 

A sandpile configuration on G is a vector r\ = (t] v ,v £F)£ z)YK The number rj v represents the 
number of grains of sand present at the vertex v in the configuration 77. When this number exceeds 
a certain threshold, the vertex is said to be unstable, and will topple, sending one grain of sand to 
each of its neighbours. Typically and throughout this paper, the threshold is set to the degree of 
that vertex. The sink plays a special role in that it can absorb any number of grains, and as such 
never topples. 

Two different configurations play an important role in this paper. The first one rjQ ax (or ry max 
when G is clear from the context) is the configuration 

r ] ^:=(d G (v),v€V), (1) 

where d G (v) is the degree of the vertex v in G. The second one is l a for some vertex a in V , 

t a :=(6 a , v ,veV), (2) 

where 5 a)V is the Kronecker symbol, meaning that all co-ordinates of t a are except for at position 
a, where it is 1. 

Definition 2.1. A configuration r] = (r] v ,v € V) is called stable if i] v < d G (v) for all v G V. We 
write Stable(G) for the set of all stable configurations on G. 

^max - g c j ear jy maximum stable configuration in terms of the total number of grains in the 
configuration. 
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Now define the toppling operator T x corresponding to a toppling at x S V by 

T X { V ) = r,-d G {x)l x + V ( 3 ) 

where configurations are added site- by-site. A toppling T x is called legal if r\ x > d (x). 

We now define for any configuration rj its stabilisation S(rj) as follows. If 77 is stable, then 
S(rj) = rj; otherwise 

S( V ) = T Xm (...(T Xl ( V ))-), (4) 

with the requirements that T Xl , ...,T Xm is a sequence of legal topplings, and that the configuration 
S(rj) is stable. The fact that this operator is well-defined, i.e. the stabilisation of a configuration 
is independent of the order of the topplings that is used to stabilise it, is not immediately obvious. 
We refer the interested reader to Proposition 3.7 (and Lemma 3.6) in Redig [T2"] . 

Markov chain structure of the ASM 

The ASM has a Markov chain structure which is defined as follows. Assume that there are defined 
some i.i.d. random variables (Xi,i > 1) taking their values in V according to a distribution ji 
(where the support of ji is V), defined on a common probability space (fl,A, V). The sandpile 
process starts from any stable configuration 770 . We define a Markov chain (rji,i > 0) with values 
in Stable(G). Given rji-i for any i > 1, rji is obtained as follows: 

• Add a grain at position Aj to the configuration rji-\. Let <q\ be the obtained configuration. 

• Let rji be the stabilisation of rj^, that is rji = S(r^) (in some cases no toppling is needed). 

Since a grain can be added anywhere with positive probability, it is immediately seen that from 
any configuration rj, the maximal configuration rj max can be reached with positive probability. It 
follows that the set of recurrent configurations for the Markov chain is the unique recurrent class 
containing the maximal configuration. We denote this set by Det(G), and as previously noted, it 
is in bijection with the set of spanning trees of G (see Redig [12J). 

2.2 The SSM 

In this model, we make the topplings random. A probability p G (0, 1] is fixed, and when a site is 
unstable, each neighbour independently has a probability p of receiving a grain from the unstable 
site. In this way, an unstable site may remain unstable after toppling, but as in the original model, 
the process continues until the configuration is stable. Different topplings are done independently. 
If p = 1, this reduces to the ASM, so we assume this is not the case. 

There are some slight difficulties in establishing that this model is indeed well-defined. Two of 



them are taken into account in the proof of Theorem 2.2 below. A third one arises by noticing that 
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topplings commute in the ASM. However, in the SSM, when two legal topplings are possible, it is 
not clear that they commute, since we only talk about what to do "in probability" . To solve this, 
we describe a probability space where the laws of the topplings are indeed those discussed at the 
beginning of this section and for which the commutation of legal topplings takes place as needed. 

We now give a formal definition of the SSM. The definition of terms introduced in relation to 
the ASM, apart from the toppling operator, are unchanged. All the random variables discussed 
below are defined on a common probability space (Q,A, V), and they are all independent. 

• The variables (A»,i > 0) are distributed according to a distribution n with support equal to 
V. As before, they represent a sequence of arrival places of grains. 

• For any x G V, e G E such that e = {x,y} and any i > 0, Bi(x,e) is a Bernoulli random 
variable with parameter p. It represents the number of particles going from x to y along e 
due to the ith toppling of x. 

We now define T x as the ith toppling at x G V 

Tf(i 1 ) = 11 -Y,B l {x,e)t x + Bi{x,e)l y . (5) 

&x y^ s ,e={x,y}eE 

We will call any realisation of such a toppling a stochastic toppling. Note that with positive 
probability, the ith toppling at v is a full toppling, meaning that Bi(x, e) = 1, Ve B x, in which case 

— (i) 

T x coincides with the ASM toppling T x . 

On the probability space (£l,A, V), the SSM is well defined when a starting configuration is 
specified. We then use the same definition of a legal toppling, and define the stabilisation of a 
configuration r/ by 

RS( V ) = T Xn (...(T xl (r l ))-), (6) 

where we take T Xj = T x . if the vertex Xj appears exactly i — 1 times in the sequence of topplings 
T Xl , ...,T Xj _ 1 , and with the same requirements as for the ASM, i.e. that T Xl , ...,T Xm is a legal 
sequence of topplings and RS(r/) is stable. The notation RS stands for the "random stabilisation" 
we define. 

Theorem 2.2. For any graph G in Q , the stabilisation operator RS is almost surely well-defined. 

Proof. Two things have to be shown. The first one is that the stabilisation process eventually ends. 
Here, as opposed to the ASM, a loop may appear in the toppling process. For example, if a ~ b 
and a and b are both unstable, and for all i we have T a = If, and T b = l a , then a repeated 
loop occurs where a only sends one grain to b and b only sends one grain to a, ad infinitum. This 
example is obviously contrived, and it is easy to see that this occurs with probability 0, but some 
more complex loops could occur. 

To show that this almost surely does not happen, we note that in the ASM, stabilisation takes a 
finite number of topplings regardless of the initial configuration. Let N v be the maximum number 
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of full topplings at v needed (by the ASM) to stabilise any configuration arising from the addition 
of 1 grain to a stable configuration; let N = max„ g yA^. Since there are a finite number of 
such configurations, N is finite. Now for any i and x, T x is the "full toppling" T x with positive 
probability. A simple renewal argument shows that in an infinite sequence of topplings, the number 
of full topplings done at every vertex is almost surely larger than N. From here the conclusion 
follows (note that this doesn't mean that full topplings are needed to stabilise a configuration). 

The second thing we have to prove is that that the stabilisation operator is well-defined; that 
is, starting from an unstable configuration r], changing the order of the site stabilisations does not 
change the final result RS(?y). This is the case on (J), A, V) since whatever is done elsewhere, the ith. 
toppling at v is T v . We can then use the same argument as in Lemma 3.6 and Proposition 3.7 of 
Redig [12], with a simple adjustment to take into account the fact that i h- > T x is not constant. □ 

Markov chain structure of the SSM 

The SSM has a Markov chain structure which is analogous to the ASM. The sandpile process starts 
from any stable configuration rjQ. We define a Markov chain (rji,i > 0) with values in Stable(G). 
Given 774—1 for any i > 1, we obtain m as follows: 

• Add a grain at position to the configuration T)i—\. Let rf i be the resulting configuration. 

• Let rji be the stabilisation of 77^, that is r\i = RS(t^) (in some cases no toppling is needed). 

Note that RS is not defined ex nihilo as in the ASM. On (Q,A, V), when RS is applied, its 
action depends on all the previous topplings taken to reach the current state. In view of this, it is 
more proper to write RS^ or even RSj(r/o) instead of RS. However, for the sake of brevity we write 
RS instead. 

Once again, since the support of /j is V, the maximal configuration rj max is recurrent, and the 
set of recurrent configurations for the Markov chain is the unique recurrent class containing rj max . 
We denote this set by Sto(G), and call any element of it stochastically recurrent (or SR). We call the 
configurations which are recurrent under the ASM deterministically recurrent or DR in contrast. 

Proposition 2.3. For any G G Q , 

Det(G) C Sto(G). 

Proof. This is easily seen because S always uses a finite number of topplings. When constructing 
a stochastic stabilisation, the deterministic stabilisation must have a positive probability of being 
reconstructed using only full topplings. □ 

The two sets Det(G) and Sto(G) are not equal in general. A simple counter-example is shown 
in Figure [TJ For this graph, the configuration which has 2 grains at each vertex is SR. It can be 
reached from the DR configuration [3, 1, 2] (with 3 grains at the vertex v\ leading to the sink), by 
adding a grain to v\ and toppling it, sending one grain to the sink and one to vi- However, it is 
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not DR as it fails the burning algorithm test — the graph with the sink removed is a forbidden 
subconfiguration. 



as 

Figure 1: A graph for which Det(G) is a strict subset of Sto(G). In this case Det(G) = 
{[3, 1, 2], [3, 2, 1], [3, 2, 2]}, and Sto(G) = {[2, 2, 2]} U Det(G). 

To determine Sto(G), it is clear that rj is SR if and only if there exists a finite sequence of 
adding of grains and topplings such that rj is reached from r\ max through this sequence. In the rest 
of the paper, we prove a more useful characterisation of the stochastically recurrent states. 

3 Main results 

In this section, we state the main results of this paper. Our two main results are a characterisation 
of the stochastically recurrent states in terms of graph orientations and a recurrence for the lacking 
polynomial (which we define below) . Proofs of these results will be given in Section |4j 

3.1 Graph orientations 

Our first result characterises the SR states in terms of graph orientations. Take a graph G = 
(V U {s}, E) 6 Q. We define an orientation on G to be an orientation of each edge of E (when G 
has multiple edges, all of them are oriented independently). We write (a, b) or a — > b to denote 
that the edge {a, 6} is oriented from a to b. 

Definition 3.1. Let G = (V U {s},E) £ Q. Take a sandpile configuration n on G. We define the 
lacking number of rj at v as the number of grains at v less than its maximum value: 



Now let O be an orientation on G and let \r\Q(v) be the number of incoming edges to v in O. We 
say that r] is compatible with O (and likewise O is compatible with rj) ifVv£ V , 





{v)-rj v . 



in>)>l + /: 



(7) 



We denote the set of stable configurations that are compatible with O as comp(O). 
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In situations where it is clear, we will omit the superscript G for brevity. 

Note that there may be several configurations compatible with a particular orientation. Like- 
wise, there may be several orientations compatible with any given configuration. For instance, the 
maximal configuration rj max is compatible with any orientation where each vertex has at least one 
incoming edge. 

Theorem 3.2. Let G = (VU{s},E) 6 Q . Then a (stable) configuration rj is stochastically recurrent 
if and only if there exists an orientation O on G such that n E comp(O). In other words, 

Sto(G) = Jcomp(O), (8) 

o 

where the union is taken over all orientations on G. 

Furthermore, there is also a way of characterising DR configurations using orientations. 

Theorem 3.3. A (stable) configuration n is deterministically recurrent if and only if there exists 
an orientation O of G with no directed cycles such that n G comp(O). 

This theorem is intuitive given the bijection between DR configurations and spanning trees, as 
any spanning tree can induce a (not necessarily unique) orientation with no directed cycles. 

3.2 The lacking polynomial 




as as 

Figure 2: The graph G on the left has Det(G) = 

{[3, 1, 2, 3], [3, 2, 1, 3], [3, 2, 2, 3], [3, 3, 1, 2], [3, 3, 1, 3], [3, 3, 2, 1], [3, 3, 2, 2], [3, 3, 2, 3]} and 
Sto(G) = Det(G) U {[1, 3, 2, 3], [2, 2, 2, 3], [2, 3, 1, 3], [2, 3, 2, 2], [2, 3, 2, 3], [3, 2, 2, 2]}. Thus the 
lacking polynomial is Lq{x) = 1 + 4x + 9x 2 . The graph on the right has lacking polynomial 
l + 4x + 8x 2 . The DR configurations are equinumerous in both graphs since they are in bijection 
with the sets of spanning trees, which can be identified with each other. This does not hold for 
the SR configurations. 



Our second result uses Theorem 3.2 in order to classify SR configurations according to the total 
number of grains, or equivalently to the number of grains removed from the maximal configuration. 
We do this by means of the lacking polynomial, which we now define. 
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Definition 3.4. Let G = (V U {s},E) G Q. The lacking polynomial L G of G is the generating 
function of the stochastically recurrent configurations on G, with x conjugate to the number of 
lacking particles in the configuration: 

L G (x)= Yl ^ ( 9 ) 

r?eSto(G) 

where 

'to) = £ 

v&V 

An example of the lacking polynomial is shown in Figure [2j Note that we can use the lacking 
polynomial to count the number of SR configurations, as |Sto(G)| = L G (1). 

Before stating the main result of this section (Theorem 3.9), which gives a Tutte-like formula 
to compute L G , we state (without proof) some propositions concerning some special graphs which 
serve as an initialisation for the computations (since L G will be expressed in terms of the lacking 
polynomials of some graphs smaller than G). Related illustrations can be found in Figure |3j 

Proposition 3.5. Let G = {V U {s}, E) G Q. 

1. IfV = {u} and there are k edges between u and s, then 

k-i 

L G (s) = £s\ 

i=0 

2. If G is a tree, then Lq(x) = 1. 

3. If we can write G as the union of connected graphs Gi = (Vi U {s}, Ei), i = 1, . . . , k, so that 
the V% are mutually disjoint, then 

k 

L G (x) = l\L Gt (x). 
i=i 

We say that G is the product of the Gi. 

We now state that the pruning of "tree branches" of a graph G does not change its lacking 
polynomial. 

Definition 3.6. Let G = {VU{s},E) G Q. A tree branch of G is a subgraph T = (V' U {r}, E') of 
G which is a tree attached to the rest ofG at the vertex r. In other words, T is a tree, d T (v) = d (y) 
for any v G V , and r G V \ V . 

Lemma 3.7. Let G = {V U {s},E) G Q, and let T = (V U {r}, E') be a tree branch of G. Define 
G\T := ((V \ V) U {s}, E \ E') as G with the tree branch removed. Then 

l g\t{%) = L G (x). 
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Figure 3: From Proposition 3.5 the first graph has lacking polynomial 1 + x + x + x , the 
second graph (which is a tree) 1, and the lacking polynomial of the third graph is the product of 
the lacking polynomials of the 3 graphs surrounded by dashed lines (all containing a sink at the 



same place). Lemma 3.7 says that in the fourth figure, the lacking polynomial is not affected by 



the removal of the two tree branches surrounded by dashed lines. 



Note that this lemma implies the second statement of Proposition 3.5 We now give a formula 
allowing one to compute the lacking polynomial for a given graph; this formula is similar to the 
one used to compute the Tutte polynomial. First we require some definitions of edge deletion 
and contraction, similar to those used in the Tutte polynomial relation (see e.g. Bernardi [1] and 
references therein). 

Definition 3.8. Let G = (V U {s},E) E Q , and consider an edge e = {x, y} £ E, with x,y ^ s. 

1. Edge deletion. The graph G\e is the graph G with e removed, i.e. G\e = (VL){s}, E\{e}). 

2. General edge contraction. Define the graph G.e as follows: 

• If e is simple, then G.e is G with e contracted, i.e. G.e = (V U {x.y, s} \ {x, y}, E \ {e}), 
where edges adjacent in G to either x or y are now connected to x.y instead. 

• If e has multiplicity k > 2, contract one of these edges as above, and replace the other 
k — 1 edges with k — 1 edges {x.y, s}. 

We illustrate these operations in Figure |4} 

Theorem 3.9. Let G = (V U {s},E) £ Q, and let e be an edge of E which is neither a bridge (i.e. 
removing e doesn't disconnect the graph), nor connected to the sink. Then 

L G (x) = xL G \ e (x) + L G , e (x). (10) 

By way of contrast, the corresponding relation for the Tutte polynomial, where e is not a loop 
or bridge, is 

Tg(x, y) = T G \ e (x, y) + T G , e (x, y), 
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Figure 4: An example of the deletion and contraction operations. G with the edge e marked is 
on the left; G\e and G.e are in the centre and right respectively. 



although here loops are allowed in G and so G.e denotes regular edge contraction rather than the 
version defined above. 



One can check that the use of Theorem 3.9, Proposition 3.5 and Lemma |3.7| allows one to 



compute Lq for any graph G in Q without referring to the SSM. More specifically, Theorem 3.9 



expresses Lq in terms of graphs with one less edge. We can continue to use this theorem, and case 3 



in Proposition 3.5, until we express Lq in terms of graphs which only contain edges to the sink and 



bridges. The bridges must then form tree branches which are removed by Lemma 3.7 and case 1 in 



Proposition 3.5 provides the lacking polynomials of the remainders. It is not immediately obvious 
that the polynomial so obtained does not depend on the edges that we choose to delete/contract; 



however, since the lacking polynomial itself is well-defined, this follows from Theorem 3.9 
A simple corollary of Theorem |3.9| gives us the degree of the lacking polynomial. 

Corollary 3.10. Let G = (V U {s}, E) 6 Q. Let the number of cycles c(G) be the number of edges 
it is necessary to remove in order to turn G into a tree, i.e. c{G) = \E\ — \ V\. Then the degree of 
the lacking polynomial Lq is the number of cycles c{G). 

Proof. We use induction on \E\. If \E\ = 1 then G = ({«, s}, {{u, s}}). Thus Lq = 1, c(G) = as 
desired. Now consider a graph G = (FU {s}, E), and an edge e which is not a bridge or connected 
to the sink. If there is no such edge, then after removing tree branches, which does not affect c(G), 
we can express G as the product of graphs of the form described in the first case of Proposition 
3.5, which obviously satisfies the corollary. Otherwise, the graphs G\e and G.e both have \E\ — 1 

c(G) — 1 and c(G.e) = c(G). The result then immediately follows by 

□ 



edges. Moreover, c(G \ e) 
induction using ( |10[ ). 
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Figure 5: After the vertex topples according to the dashed arrows, only the edge e is reversed. 

4 Proofs 

4.1 Proof of Theorem O 

Let 

V(G) = Jcomp(O) 
o 

be the set of stable configurations on G compatible with some orientation. We wish to show that 
Sto(G) = V(G). We do this by showing that they are subsets of each other. 

Lemma 4.1. Let G = (V U {s}, E) G Q. We have 

Sto(G) C V(G). 

Proof. Firstly we show that rj max G V(G). For this, consider a spanning tree of G and orient all 
edges of this tree outwards from the sink; that is, if {a, b} is an edge of the spanning tree such that 
a is closer to the sink than b, orient the edge a — > b. Orient all remaining edges in any direction. 
Then for the resulting orientation every non-sink vertex has at least one incoming edge. Condition 
([7]) now shows that rj max is compatible with this orientation, so rj max G V(G). 

Now let O be an orientation compatible with r] max , and let rj G Sto(G). We shall construct an 
orientation O' which is compatible with rf . To do this, consider a history of grain additions and 
stochastic topplings that leads from rj max to rf. We construct O' iteratively from O by making the 
following changes to the orientation at each step of this history: 

• If a grain topples from a vertex a to a neighbour b and the edge {a, b} is oriented a — > b, we 
reverse the orientation of this edge so that it is oriented b —> a. 

• Otherwise — that is, if the edge {a, 6} is oriented b — > a, or if a grain is added — do nothing. 

These changes are shown in Figure |4.1[ 

We now show that rf G comp(0 / ). Let 771 = r] max ,r]2, ■ ■ ■ ,r]k = if be the configurations con- 
structed at each step of the history, some of which may be unstable. Let 0\ = O, Ok = O' be the 
corresponding orientations. We show by induction on i that for any i G {1, . . . , k}, rji is compatible 
with Oj. For i = 1 we have rj max G comp(O) by definition. At any fixed i > 1, there are two cases: 
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1. r\i is reached from T)i—\ by addition of a grain at some vertex a. 

For any vertex v, we have l Vi (v) < l rH _ 1 (v) and by construction \r\Oi(v) = inOj_i(^)- Since 
is compatible with Oj_i, condition ([7]) is satisfied for rji and Oj for all vertices in V. 
Thus rji is compatible with Oj. 

2. 7y, is reached from rji-i through a (legal) toppling of some vertex a. 

Let the neighbours of a be bi, . . . , b m , and suppose the toppling transfers a\, . . . , a m grains 
from a to each of these vertices respectively. Since the toppling is legal, ryj-i| a > d(a) + 1, and 
so after the toppling, l m (a) < . — But now grains have toppled out from a along ^ ■ dr- 
edges, so by construction all these edges are oriented towards a. Therefore ino^a) > X)j a j 
and condition ([7]) is satisfied at a. 

It remains to check that (JjJ is still satisfied at each bj. At most aj incoming edges to bj in O l -\ 
are now outgoing in Oi, so ino^j) > ino^^j) — «j- Furthermore, l Vi (bj) = — ctj, 

so by induction ([JJ is satisfied at bj. 

No other vertices apart from a and its neighbours are changed, so rji is compatible with Oi. 

Taking i = k, we have shown that O' is compatible with r/ . This completes the proof of the 
lemma. □ 

Lemma 4.2. Let G = (V U {s}, £) G £. We have 

V(G) C Sto(G). 

Proof. We use induction on the number of vertices |V|. If \V\ = 1, then G is simply one vertex 
v , connected to s by say k edges. Then V(G) consists of all configurations r] where l-qiv) < k — 1, 



which are SR by Proposition 3.5 



Suppose now that the lemma is true for any graph with \V\ = k — 1. Take a graph G = 
(V U {s},E) with \V\ = k and a configuration 77 G V(G). We wish to show that ry is SR. To do 
this, take an orientation O of G such that i] G comp(O), and choose a vertex n connected to s. We 
can assume that all edges connecting u to s are oriented towards m (if some are not, we can reverse 
them and the resulting orientation remains compatible with 77). Write m u for the number of such 
edges. 

Let G u be the graph obtained from G by identifying u with s, removing all {u, s} edges (denote 
the new sink by s.u). Note in particular that the degree of any vertex in G u (except s.u) is equal 
to its degree in G. We now let O u be the orientation on G u coinciding with O on all edges still 
present, except that any edge connected to s.u is oriented away from it. 

For any vertex v in G u , let k v be the number of edges oriented v — > u in O (k v > 0). We 
define a configuration r\ u on G u by rj% = rj v — k v for v 7^ s.u in G u . Now for any such v, we have 
in^(w) = inQ(v) + k v and l^u(v) = l^(v) + k v . Since r/ G comp(O), we deduce that rj u G comp(0 M ). 
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Hence by induction r] u E Sto(G u ), and therefore there exists a history of grain additions and legal 
topplings that leads from r] max to r\ u on G u . 

We now start from rj max on G and copy this history (all legal topplings remain legal because 
the degrees in the two graphs are identical). This results in a configuration rj on G with rj' v = rj% 
if v ^ u. We then continue the history by either adding grains to u or repeatedly toppling grains 
from u to s until u has d ■ (it) + 1 grains, i.e. it is minimally unstable. 

We now make one final toppling at u, sending k v grains to each of its neighbours v, and 
max(^f(u) — k v + 1,0) grains to the sink. In order to do this, we must have at least l9{u) — 
J2 V k v +l edges {u, s}. However, since rj and O are compatible, we know that l~*(u)+l < m u +E v k v , 
so this is true. Denote by rj' the configuration we finally reach. 

We have: 

• If v is neither u nor one of its neighbours, then Gi{v) = l3(v) = l^u{v) = l^(v). 

• If v is a neighbour of u, then l^„(v) = l3(v) — k v = l^u(v) — k v = l^(v). 

• >Ev k v + - Ev k v + 1 - 1 = 

Together, this shows that 77" is identical to r\ except at u, where it may have less grains. We then 
merely add the difference in grains to u, and have thus created a history of grain additions and 
legal topplings which leads from r\ max to 77. Therefore r/ is SR and the lemma, and Theorem 3.2 is 



proved. □ 
4.2 Proof of Theorem B 

Let rj € Stable(G) be compatible with an orientation with no directed cycles. Let O be such an 
orientation. Since O has no directed cycles, and we can take all edges adjacent to the sink as 
oriented away from it, we can order the vertices of G as s = Vq, vi, . . . , v n so that there exist no 
edges Vj — > V{ for i < j. 

Now we apply the burning algorithm to 77. We claim by induction that this can burn the vertices 
in the order described above. Since s = vq, the initial condition is trivial. Now suppose we have 
burned vertices Vq, . . . All incoming edges to V{ are burnt, so the number of unburned edges 

adjacent to Vi is d(vi) — \r\o(vi). But by ([7]), 

rj Vi = d(vi) - l v (vi) > d(vi) - ino(v») + 1. 

Therefore Vi can be burnt. Thus the burning algorithm burns all the vertices of the graph, and 77 
is deterministically recurrent. 

Conversely, let 77 € Det(G). Now apply the burning algorithm to 77, and every time we burn 
a vertex v, orient all edges from previously burnt vertices to v as incoming edges to v. Since all 
vertices are burnt, this produces a full orientation O on G, which obviously has no directed cycles. 
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From the burning condition, we know that for all v, rj v is greater than the number of unburnt edges, 
which is d(v) — \x\o{v). This gives 



Since these are integers, this means that ([7]) is fulfilled for all vertices. Thus O is compatible with 
r\ and the theorem is proved. □ 

4.3 Proof of Lemma 13.71 

We construct a bijection <3> from Sto(G\T) to Sto(G) such that for any ry £ Sto(G\T), l(rj) = l(®(r])). 
Define for rj £ Sto(G \ T) and v £ V 



Take an orientation O on G \ T which is compatible with r/, and extend this orientation to G 
by orienting each edge in T away from r. Then is compatible with the resulting orientation, 
so <&(r]) £ Sto(G). Moreover, <!> is clearly an injection. 

It remains to show that it is surjective. To see this, consider a configuration rj £ Sto(G) and 
a compatible orientation O. Each vertex in V must have at least one incoming edge in E' . But 
since T is a tree, \V'\ + 1 = \E'\ + 1. Thus each edge in E 1 points to a different vertex in T, so 
' n o( v ) = 1 f° r an v ^ V • This implies that l~(v) = 0. Furthermore, all edges in E' adjacent to r 
point away from it. 

Now define rj on G\T according to fj^ 1 \v) = l^j{v), and let O' be the restriction of O to G\T. 
We have \n ( ^' T (v) = in^w) for all v € G \ T, so clearly rj is compatible with O' and $(?/) = r\. 
Thus $ is a bijection. 

4.4 Proof of Theorem 13.91 

Let G = (VU {s}, E) S Q, and let e be an edge of E which is neither a bridge nor connected to the 
sink. Write e = {a, b}. For r/ € Sto(G) we distinguish the following two cases: 

(A) There exists an orientation O on G, compatible with rj, such that e is oriented a — > b in O, 
and l®{b) > 0. 

(B) For all orientations O compatible with rj, all {a, b} edges are oriented b — > a, or l^(b) = 0. 
We write rj £ StO( x )(G) if rj satisfies condition x £ {A,B}. Obviously Sto(G) = StO(^)(G) U 



l v (v) = d(v) - 7) v < d(v) - (d(v) - ino(w)) = ino(v)- 




if v £ V, 
l^ T (v) otherwise. 



Sto (B) (G). 

Now we define a function / : Sto(G) — > Stable(G \ e) U Stable(G.e) as follows: 



• If rj £ Sto^^G) then f(rj) = r/ — t a is a configuration on G \ e. 
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• If 77 G Storm (G) then f{rf)\ v = J ^ a ^ b ' ' is a configuration on G.e. 

I 7] v otherwise, 

To simplify further calculations, we note that in the first case, ^^\( a ) = ^(°) an d ^fM^) = 
l%(b) - 1. In the second case, lf;\(a.b) = l^{a) + l%(b). 

It is easy to see that if n is stable, f(n) is also stable — the lacking number can decrease by at 
most one, and this occurs only at b when rj G StO(^)(G), where by definition l9(b) > 0. Moreover, 
we have 

«*(„)) = ! M- 1 if ^Sto (A) (G), 
1 Kv) ifr/GSto (B) (G). 



In light of this, it is sufficient to show the following theorem to prove Theorem 3.9 since the two 
sets Sto(G \ e) and Sto(G.e) are disjoint, being configurations on different graphs. 

Theorem 4.3. Let G = (VU {s}, E) 6 Q . The function f defined above is a bijection from Sto(G) 
to Sto(G\e)uSto(G.e). 

This theorem is itself a direct consequence of the four following lemmas. The first lemma shows 
that for any n £ Sto(G), the configuration f(jj) is indeed stochastically recurrent. 

Lemma 4.4. For any G = (V Li {s}, E) in Q, 

/(Sto(G)) C Sto(G \ e) U Sto(G.e). 

Proof. Take r\ £ Sto(G), and fix an orientation O on G such that n G comp(O). There are two 
cases. 

1. 7 ? GSto (A) (G). 

By construction, we may choose O such that e is oriented a — > b. Let O' be the orientation 
on G \ e which is identical to O on all edges of G \ e. Then for any vertex v, 

ing}» = ing(v) - 5 Vyb > 1 + l°{v) - 8 v>h = 1 + l^(v), 



so f(rj) G comp(0 / )- Thus, by Theorem [372| /(ry) G Sto(G \ e). 
2. n G Sto (B) (G). 

Let k be the number of {a, b} edges in G. Then in G.e, these will be replaced by k — 1 edges 
{a. b, s}. Orient these as s — > a.b, and orient all other edges of G.e as they are oriented in O. 
Denote by O' the resulting orientation on G.e. Then 

ing; e (a.6) = ing(a) + ing(6) - 1 > (Z«(o) + 1) + (l%(b) + 1) - 1 = lf(^(a.b) + 1, 

so condition (M) is satisfied at a.b. Since ing; e (t;) = ing(u) for v ^ a.b, it is clearly also satisfied 



elsewhere, so f(n) G comp(0 / ), and by Theorem 3.2, f(rj) G Sto(G.e). 



16 




lb 




lb 



Figure 6: Counting the number of edges in O' (left) and O (right). The {a, b} edges are 
oriented according to the value of i a and l^ e {a.b). 



□ 

We write /m) (resp. fm)) f° r the restriction of / to the set StO(^)(G) (resp. Stcvm(G)). We 
will show that each of these are bijections onto their respective images. 

Lemma 4.5. For any G = (V U {s},E) in Q, the function fr^\ is a bijection from Stcv^G) to 
Sto(G\e). 

Proof. The fact that f^ is injective follows immediately from the definition of /. To show that 
/m) is surjective onto Sto(G \ e), let rj G Sto(G \ e) and take an orientation O on G \ e compatible 
with v. Let O' = O U {a — > b}, and define n' = n + t a as a configuration on C It is obvious that 
n' is stable and f(rj') = n. Then for any vertex v, 



ing,(u) 



in 



GV{ v) + 5 v , b >l + l^ e (v) + 5. 



o 



v,b 



Thus r/ G comp(O') so n' is SR. Moreover, e is oriented a — > b in O' and 6(6) > 0, so 77' € Sto (A) (G) 



as desired. 



□ 



Lemma 4.6. For any G = (V U {s},.E) in £/, the function f^) is surjective onto Sto(G.e). 

Proof. For any rj G Sto(G.e), let O be a compatible orientation on G.e. Let be the multiplicity 
of the edge e in G (as before we may have k = 1). We may assume that the k — 1 edges {s, a. 6} 
in G.e are all oriented s —> a.b in O. Write i a (resp. %) for the number of edges oriented into a. b 
in O which correspond to edges into a (resp. b) in G, from vertices other than b (resp. a). This is 
illustrated in Figure [6| 

Since n is compatible with O, we have 



in 



G.e, 
O 1 



a.b) = i a + ib + k-l>l + l^- e (a.b). 



(11) 



Let /(b) : Sto(G) — > Stable(G.e) be the extension of f^B) to all of Sto(G), i.e. 

f(B){v)\v = 



Va + rib- 2 if v = a.b, 
rj v otherwise. 
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The difference between fm) an d f(B) lies i n the fact that /(#) is defined over all Sto(G), whereas f(B) 
is defined only on StO(g)(G). We will now define a configuration rf G Sto(G) such that f(B)(v') = 7 7- 
Firstly let rj v = r) v if u ^ a, 6, a. 6. Likewise let O' be an orientation on G where all edges in G.e are 
oriented identically to O (the remaining edges are as yet unspecified) . Obviously Q is satisfied for 
7/ and O' at vertices other than a and 6. We now assign grains to a and b, and orientations to the 
{a, b} edges, according to 3 cases. 

1. i a = 0. 

We set l^,(a) = 0, = l^' e {a.b). In O', orient one {a, b} edge as 6 — > a and the remaining 

as a —7- 6. Clearly, ([7]) is satisfied at a (which has one incoming edge), and 

ing,(6) = i b + /c-l = i a + i fc + A;-l>l + lf e {a.b) = 1 + §(6), 

from (11), so ([7]) is also satisfied at b. 



2. 1 < i a < l + lf e (a.b). 

We set l%(a) = i a - 1 > 0, §(6) = l^ e (a.b) - l%(a) > 0. In O', orient all A; edges {a, b} as 
a — > b. Then 

ing,(a) = i a = l + l%(a), 
ing,(6) = i 6 + fc > 1 + tf- e (a.b) - i a + 1 = 1 + §(&), 
again using (11), so ([7]) is satisfied at both a and 6. 



3. i > l + lf e {a.b). 

We set Z^(a) = l^ e (a.b), l^,{b) = 0. In O', orient all k edges {a, b} as a — > b. Then we have 
ing/(a) = i a > 1 + /^(a), and ing/(6) = % + k > 1 = 1 + l^,(b). Thus condition (jT]) is satisfied 
at a and b. 

Now, in each of these cases we have l^ e (a.b) = l^,(a) + l?r(b), so f(B)iv') = Likewise, 77' is 
compatible with O' , so 77' € Sto(G). It remains to show that we may choose an 7/ G StO( B )(G) so 
that /(t/) = 77. 

To show this, define V b nm = min{/y(6) : r( G Sto(G) s.t. /( j b)( 7 ?') = ^l- Since there exists at 
least one such 7;', this is well-defined. Now take 7/ G Sto(G) such that /^(fr) = Z™ n and /(b)(V) = r 7- 
We show that 77' G StO( B )(G). 

If Z™ n = 0, this is true by definition. Now suppose that I7 im > and there exists an orientation 
O' on G compatible with 77' with an edge oriented a — > b. We define the orientation O" as O' 
with that edge reversed and all other edges oriented as in O. Likewise, define the configuration 
77" =r f + l b - t a , so that §,(0) = §(a) + 1, l^„(b) = l^(b) - 1 > 0, and l%,(v) = l^(v) elsewhere. 
Since 77' G comp(O'), we have 77" G comp(0") by construction. Now we have /(_b)( 7 ?") = V an d 
77" G Sto(G), but l^//(b) < l™ m . This is a contradiction of the definition of Z™ n . Therefore no such 
orientation O' exists, and 77' G StO(£)(G). □ 
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Lemma 4.7. For any G = (V U {s},E) in Q, the function /Vm ^ s injective. 

Proof. Let 771,7/2 G Sto (B) (G) such that /(771) = f( m ), that is Zg(a) + Z{*(6) = ig(a) + Zg(6) 
and = ^(v) if f / a, 6, but suppose 771 7^ 772- Assume without loss of generality that 

l„ 2 (b) > iS(6) > 0. Now choose compatible orientations Oi and O2 respectively on G. Suppose 
that the edge e has multiplicity k > 1 in G. Since 772 G StO(#)(G) and > 0, these k edges 

must be oriented b — > a in O2 ■ 

Now, if ing 2 (a) > 1 + lg 2 (a), then reversing the orientation of e to a — > b in O2 results in 
another orientation compatible with 772, contradicting the fact that 7/2 £ StO(m(G). Therefore 
these quantities are equal and 

ing^a) > 1 + Ig (a) > 1 + lg(a) = ing 2 (a). (12) 

Now let A be the set of edges in O2 which are oriented differently from 0\. We define a subgraph 
T = (V(T), E(T)) of G as the union of all directed paths in O2 starting from a whose edges are in 



A (and the induced vertices). From (12), this contains at least one edge adjacent to a. 

Firstly, we claim that b £ V(T). Otherwise, there exists a directed path from a to b in 02- 
Starting from O2, we may reverse e and all edges of this path to reach an orientation with the 
same number of incoming edges at each vertex as O2, and therefore compatible with 772, but with 
e oriented a — > b. This contradicts the assumption that 772 G Storm (GQ. 

Now start from O2 and reverse the orientation of e and all edges in E(T). Denote this orientation 
by (9 2 . We show that 0' 2 is compatible with 772 by checking condition ([7]) at b and vertices in V(T) 
(which include a): 

• ing, (b) = ing 2 (6) + 1 > 1 + l^ 2 (b), since e has been reversed and b £ V(T). 

• For v G V(T), all incoming edges to v in 0\ are identically oriented in 0' 2 by construction. 
Therefore ing^(v) > ing 1 (u) > 1 + lg{v) > 1 + lg 2 (v), where the last inequality is strict if 
v = a and an equality otherwise. 

This gives us an orientation compatible with 772 containing an edge a — > b. Again, this is a 
contradiction of the assumption that 772 G StO(m(G). Thus there cannot exist configurations 
Vi 7^ in StO(m(G) such that /(771) = f{rj2), and /(b) is injective. □ 



5 Conclusion 

In this paper, we have devised a generalisation of the ASM in which the topplings are stochastic. 
This model behaves qualitatively differently to the established ASM of Dhar. In particular, the set of 
recurrent states of this model contains that of the former model. We have proved a characterisation 
of these states using graph orientations. We also define a generating function of these states which 
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counts the number of "lacking" grains, and show that this "lacking polynomial" satisfies a recurrence 
relation which resembles that of the Tutte polynomial. 

There are two directions in which to advance this work. Given the many combinatorial interpre- 
tations of the Tutte polynomial, it would be of interest to see if the lacking polynomial demonstrates 
similar interpretations. In other words, the lacking polynomial may count certain combinatorial 
objects for given values of its parameter, and we would like to determine what these objects are. 
Alternatively, the lacking polynomial may be related in some way to the Tutte polynomial, and the 
nature of this relation should be determined precisely. 

The other topic of interest is to probe further into the behaviour of the Markov chain structure, 
more specifically the steady state. In the ASM, all recurrent states are equally likely, but this is 
not the case for the SSM. It would be interesting to calculate the probabilities for the stochastically 
recurrent states. Once we have done so, we can analyse the behaviour of the model in the steady 
state, and see if it displays a similar power-law behaviour to that observed for the classic model. 
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